Gaseous NO2 induces various envelope alterations in Pseudomonas fluorescens MFAF76a

Anthropogenic atmospheric pollution and immune response regularly expose bacteria to toxic nitrogen oxides such as NO• and NO2. These reactive molecules can damage a wide variety of biomolecules such as DNA, proteins and lipids. Several components of the bacterial envelope are susceptible to be damaged by reactive nitrogen species. Furthermore, the hydrophobic core of the membranes favors the reactivity of nitrogen oxides with other molecules, making membranes an important factor in the chemistry of nitrosative stress. Since bacteria are often exposed to endogenous or exogenous nitrogen oxides, they have acquired protection mechanisms against the deleterious effects of these molecules. By exposing bacteria to gaseous NO2, this work aims to analyze the physiological effects of NO2 on the cell envelope of the airborne bacterium Pseudomonas fluorescens MFAF76a and its potential adaptive responses. Electron microscopy showed that exposure to NO2 leads to morphological alterations of the cell envelope. Furthermore, the proteomic profiling data revealed that these cell envelope alterations might be partly explained by modifications of the synthesis pathways of multiple cell envelope components, such as peptidoglycan, lipid A, and phospholipids. Together these results provide important insights into the potential adaptive responses to NO2 exposure in P. fluorescens MFAF76a needing further investigations.

www.nature.com/scientificreports/ In the environment, microorganisms can be exposed to RNS from various sources, such as the atmosphere or the immune response, as well as collateral products of their own metabolism 11,12 . Some bacterial species possess denitrification pathways allowing them to use these species as terminal electron acceptors for their energetic metabolism [12][13][14] . The mechanisms developed by bacteria to counteract the effects of RNS are tightly linked to their virulence. Indeed, during the primary immune response, macrophages produce a burst of RNS and reactive oxygen species (ROS) to eliminate pathogens 11 . Thus, resistance to RNS is often a crucial part of the infectious process 15,16 .The first bacterial component to encounter exogenous RNS is the bacterial envelope. In Gram negative bacteria, this envelope is constituted by two lipid bilayers separated by the periplasmic space containing the peptidoglycan layer. The stability and functions of these membranes are the result of complex interactions between lipids and proteins. Indeed, membrane proteins functionality is often dependent on the steric hindrance caused by nearby lipids 17,18 . Membranes play an important role in the context of nitrosative stress 4 . NO • and NO 2 , the most common RNS found in the environment, is more soluble in the hydrophobic core of the membranes than in water 4 . Their great concentration in the membrane enhances their reactivity with other RNS and reactive oxygen species (ROS) such as O 2 •− , which reacts with NO • to form the extremely reactive ONOO −19 . Furthermore, NO 2 can react with fatty acids possessing at least two unsaturations in several ways such as peroxidation, nitration and nitrosylation 4 .
Pseudomonas fluorescens is considered a ubiquitous bacterium with a very adaptable metabolism. The P. fluorescens MFAF76a strain was isolated from air and possesses several virulence factors, such as the secretion of lipases and proteases, the ability to form biofilms at 37 °C and a high lytic activity on the human pulmonary cell line A549 20 . These characteristics make the strain MFAF76a a good model to study the effects of atmospheric NO 2 pollution on bacteria. The objectives of this study were to assess the physiological impact of gaseous NO 2 on the bacterial cell envelope and to characterize the potential molecular mechanisms underlying the adaptive responses to NO 2 exposure in P. fluorescens MFAF76a.

Materials and methods
Bacterial strain. The bacterial strain used in this study is the Gram negative P. fluorescens MFAF76a, an airborne strain isolated from dust clouds of the Rouen harbor installations (Normandy, France) generated during crop ship loading 20 . Exposition to NO 2 . P. fluorescens strain MFAF76a was precultured at 28 °C for 24 h in Lysogeny Broth (LB) (yeast extract 5 g/L, peptone 10 g/L, NaCl 5 g/L). Precultures were used to inoculate strain MFAF76a at a final optical density at 580 nm (OD 580 ) of 0.08 in Davis Medium Broth (DMB) supplemented with 2.16 g/L of glucose 21 . Culture was incubated at 28 °C during 16 h with agitation at 180 rpm. Bacteria were then centrifuged at 13,000 g for 10 min at 4 °C and resuspended in DMB at an OD 580 nm value of 4. Then, 125 µL of bacteria were deposited on nitrocellulose filters (Pore size 0.2 µm, diameter 47 mm, Sartorius Stedim Biotech GmbH, Germany) themselves deposited on Petri dishes containing DMB and 15 g/L agar. Bacteria were incubated 4 h at 28 °C and the filters were deposited on 15 g/mL agar plates, and placed into the exposition system described by Kondakova and coworkers 21 . Bacteria were exposed to 45 ppm of gaseous NO 2 at a N 2 /O 2 ratio of 8/2 (v/v) or synthetic air with the same N 2 /O 2 ratio for 2 h with a constant gas stream of 2 L/min. Exposition to NO 2 and air were limited to 2 h to prevent excessive drying of the cells, as well as for time constraints. After exposure, bacteria were resuspended in sterile physiological water (NaCl 9 g/L). The collected exposed (45 ppm NO 2 ) and control (synthetic air) bacterial suspensions were then used for all subsequent experiments.
Electronic microscopy. The sample preparation method for electronic microscopy was adapted from Bergeau et al. 22 . After exposition to NO 2 or synthetic air, 1 mL of bacterial suspension was centrifuged 5 min at 8,000 g at room temperature (RT). Supernatant was discarded and bacteria were resuspended into 100 µL of PBS. 50 µL were deposited onto electron microscopy grids for 5 min. The suspension was discarded, and the cells were fixated for 20 min with 50 µL of fixation solution (Paraformaldehyde 2%, Glutaraldehyde 0.2%, PBS). Bacteria were then washed three times by PBS for 5 min and washed three times with deionized water. Fixated bacteria were then observed by scanning electronic microscopy (JEOL JCM-6000, Tokyo, Japan). Bacteria with apparent envelope alterations were manually counted using ImageJ software version 1.52a (5 images were taken per experiment, N = 3), and their number was divided by the total number of cells 23 . Only cells that were fully visible on the pictures were taken into account.
Flow cytometry assay. The exposed and control bacterial suspensions after exposure were stained for 15 min with propidium iodine 30 nM (PI) and SYTO9™ 5.01 nM using the Live/Dead BacLight kit (L-7012, Thermofisher, United States) according to the manufacturer instructions. Bacteria were observed using flow cytometry (CytoFLEX S, Beckman Coulter, United States) and the Cytexpert v1.2 software. PI and SYTO9™ were excited at a wavelength of 488 nm and their fluorescence emissions were detected at 690 ± 50 nm and 525 ± 40 nm, respectively. A total of 10,000 events were collected for each measurement at a low flow (10 µL/ min). Positive control for membrane integrity was performed on bacteria by treating bacteria with 50% ethanol for 10 min prior to staining. The software was calibrated using the positive integrity control.
Whole proteome analysis. The bacterial suspension after exposition was centrifuged at 13,000 g at 4 °C for 20 min. Pellets were then resuspended in 20 mM Tris-HCl buffer (pH 7.4) and sonicated. Bacterial lysate was centrifuged at 10,000 g at 4 °C for 10 min, and their protein concentration was determined by Bradford test (Biorad™, United States). Total proteins were identified as described previously 24 . Briefly, proteins were loaded onto an SDS-PAGE stacking gel (7%). After short electrophoresis migration, proteins were reduced, alkylated . Samples were injected onto an enrichment column (C18 Acclaim PepMap100, Thermo Scientific, United States) and eluted using a three-step linear gradient (from 2 to 40% B over 75 min, from 40 to 80% B in 4 min, and 11 min at 80% B) at a flow rate of 300 nL/min. The samples were analyzed using CID (collision induced dissociation) method. All measurements in the Orbitrap analyzer were performed with on-the-fly internal recalibration (lock mass) at m/z 445.12002 (polydimethylcyclosiloxane). After MS analysis, raw data were imported in Progenesis LC-MS software (Nonlinear Dynamics, Netherlands). For comparison, one sample was set as a reference and the retention times of all other samples within the experiment were aligned. After the statistical analysis was performed, peptide features presenting a p-value and a q-value less than 0.05, and a power greater than 0.8 were retained. MS/MS spectra from selected peptides were exported for peptide identification with Mascot (Matrix Science, United Kingdom) against the database restricted to P. fluorescens A506 (http:// www. pseud omonas. com).
Membrane fluidity assay. The exposed (NO 2  where I represents fluorescence intensity, and "v" and "h" respectively represent vertical and horizontal settings of the excitation and emission polarizers. The instrumental correction factor G was calculated separately for each experiment by Eq. (2):

Quantitative reverse-transcription real time polymerase chain reaction (RT-qPCR). Total
RNAs were isolated either from pellets of bacteria exposed (45 ppm NO 2 ) or control (synthetic air) using the hot acid-phenol method 25 . Residual genomic DNA was eliminated using Turbo DNA-free kit (Invitrogen, United States). RNAs were non-specifically converted to single stranded cDNAs using the High-Capacity cDNA Reverse Transcription Kit following manufacturer's protocol (Applied Biosystems, United States). RT-qPCR experiments were performed as previously described 26 using the primers listed in Table 1. The relative quantification of the mRNAs was obtained by the comparative CT (2 -ΔΔCT ) method as described in 27 . Measurement of recA gene expression was used as an internal standard. recA has stable expression within the samples to be compared, regardless of physiological or experimental conditions 28 .  GAA GGC TGG TAC TGG TT  GAG AAT ATC GAA GAA GCG GAACA   pgsA  MT815581  GCA GGT GGC CGT TTC TAA TC  AAT CAC CAG CGC GAG CAT   clsA  MT815582  ACG GCA CCG CGA CCTT  TAC GGC CTC CTT CGC TTG T   clsB  MT815583  TGA GCC TGT CGC TGA ATC TG  AGT GGT TCT GGC TGA GGT CTTC   cdsA  MT815584  AAG TTC GGC AAG CGC AAA  GTA AAC GCC TTC CCA  www.nature.com/scientificreports/ Statistical analyses. All experiments were repeated at least three times. Mann-Whitney test was used to determine the significance between mean values of the cell sizes obtained by flow cytometry. Two tailed unpaired t-test was used to determine the significances of differences between mean values for the live/dead cytometry, fluorescence anisotropy and RT-qPCR experiments. Significance was set at p < 0.05 (*), p < 0.01(**) and p < 0.001 (***). For the proteomic study, after alignment and normalization, statistical analysis was performed using oneway analysis of variance (ANOVA) calculations. The significance of the peroxynitrite labelling was determined using two-way ANOVA calculations.

NO 2 exposure alters P. fluorescens MFAF76a cells morphology. Previous work by T. Kondakova
and coworkers determined that P. fluorescens MFAF76a showed significant modifications when exposed to 45 ppm of NO 2 , and we therefore chose to focus on this concentration 21 . To determine if NO 2 exposure changes the cell morphology, P. fluorescens strain MFAF76a cells were observed by scanning electron microscopy after exposure to 45 ppm of NO 2 (Fig. 1a). Around 35% of bacteria exposed to NO 2 displayed impaired cell surface as compared to the control (Fig. 1a). The aspect of altered cell envelopes was rough and bumpy, presenting asperities and potential cavities. Flow cytometry assay analyses showed that P. fluorescens MFAF76a cells exposed to NO 2 had a significant increase of the forward scatter (FSC), which indicates an augmentation of their size (Fig. 1b,c).  (Fig. 2b). The proteomic study data revealed variations in the quantity of proteins involved in the synthesis of uridine diphosphate N-acetylglucosamine (UDP-GlcNAc), the precursor of peptidoglycan and lipopolysaccharides (Fig. 2a). These proteins are GlmS, GlmM and GlmU. The results also revealed decreased amounts of FtsI, DacA and MurC, all involved in peptidoglycan synthesis and cross-linking after exposure to NO 2 .
The proteome study analysis also revealed important amount modifications variations of the proteins MraZ, FtsZ, FtsE, FtsA, FtsI and MinE playing a role in the formation of the septal ring 29 . NO 2 exposure alters P. fluorescens MFAF76a membrane fluidity and integrity. Aside from peptidoglycan, the other critical components in the structure of the cell envelope are the cell membranes. To clarify whether the effects of NO 2 observed by electronic microscopy and flow cytometry are also mediated by membrane alterations, membrane integrity and membrane fluidity were assessed for bacteria exposed to NO 2 or synthetic air. Membrane integrity assays were performed using the green-fluorescent membrane-permeant DNA probe SYTO9™ and the red-fluorescent membrane impermeant DNA probe propidium iodine (PI). Their fluo- www.nature.com/scientificreports/ rescence was measured by flow cytometry (Fig. 3b). Three populations of cells were distinguishable after labelling (Fig. 3a). Cells with intact membranes were only labeled by SYTO9 ™, while cells with a damaged membrane were labelled by both PI and SYTO9™. However, cells labelled by PI did not form a single population. Some cells formed a highly homogenous population displaying high red fluorescence and a green fluorescence in the lower ranges of SYTO9™ labelled cells (damaged membranes, Fig. 3a). This profile was similar to the control cells permeabilized by 50% ethanol, suggesting that their membranes completely lost their integrity. However, another, less homogenous, population of cells was fluorescent in the red and displayed a higher and more heterogeneous green fluorescence (partially damaged membranes, Fig. 3a). This profile was not found in the permeabilized cell control, suggesting that the membranes of these cells were only partially damaged. The PI fluorescence threshold used to discriminate between the partially damaged and the intact cells was put at the maximum PI fluorescence detected in non-labelled cells, sometimes resulting in the splitting of this population between partially damaged membranes and intact membranes. The SYTO9™ fluorescence threshold to discriminate between partially and totally damaged membranes was determined thanks to the cells permeabilized by 50% ethanol. Interestingly, the results showed a dramatic increase of the number of cells labelled by PI upon exposure to NO 2 . Indeed, 42.5% of the cells exposed to NO 2 were damaged, as compared to 2.5% of the control cells (Fig. 3b). Meanwhile, 18.2% of the exposed bacteria and 7.1% of the control bacteria (Fig. 3b) displayed the profile associated with partially damaged membranes. Membrane fluidity of P. fluorescens MFAF76a was determined by measuring the fluorescence anisotropy of the 1,6-Diphenyl-1,3,5-hexatriene (DPH) fluorescent probe (Fig. 3c). DPH is hydrophobic and localizes at the core of the phospholipid bilayer, at the level of the fatty acid chains 30 . Excitation of the probe by polarized light quantifies the polarization of emitted fluorescence. The polarization of emitted fluorescence www.nature.com/scientificreports/ depends on the capacity of the probe to rotate upon itself. Therefore, the fluorescence polarization is proportional to the fluidity of the lipid bilayer in which the probe is inserted 31,32 . DPH behavior can also be affected by its ratio to lipids, but previous gas chromatography analysis did not show any difference in terms of quantity, composition and cis/trans ratio of fatty acid between the two conditions 33 . Our results showed that fluorescence anisotropy was significantly reduced to 0.192 after bacterial exposure to NO 2 (45 ppm) as compared to 0.217 for control condition (synthetic air), indicating an increased membrane fluidity (Fig. 3c).
To clarify the consequences of the membrane alterations observed by cytometry and anisotropy, a proteomic profiling analysis of the proteins involved in the synthesis of membrane components was performed (Fig. 4). KdsC and LpxC, two proteins involved in the early steps of lipopolysaccharides biosynthesis were more detected after exposure to NO 2 . Furthermore, three proteins involved in the formation of the acetyl-CoA-carboxylase complex (ACC), responsible for the early step of fatty acid biosynthesis, were present in lower amounts in the condition exposed to NO 2 . These proteins are AccA, AccC and AccD. Furthermore, FabG and FabA, involved in latter steps of fatty acid biosynthesis were overproduced after NO 2 exposure. Finally, FadE, involved in the first step of fatty acid metabolism through beta-oxidation was found in lower amounts compared to the control condition. However, proteins involved in phospholipid biosynthesis were not identified.
In a previous study, we showed slight modifications in the phospholipid composition of P. fluorescens MFAF76a membranes after exposure to NO 2 , including the disappearance of a minor unknown glycerophospholipid (UGP), but no lipid peroxidation seems to take place 33 . These previous results combined with the observed membrane alterations noted in Fig. 1 suggest that NO 2 exposure could have an impact on the cell phospholipids. By searching the genome of MFAF76a strain, ten genes involved in the different steps of phospholipid biosynthesis from glycerol-3-phosphate (glycerol-3P) and acetyl-coenzyme A (acetyl-CoA) were found (cdsA, clsA, clsB, clsC, pgpA, pgsA, pssA, plsC, psd and pcs) (Fig. S1a). Since the proteins involved in phospholipid www.nature.com/scientificreports/ biosynthesis in P. fluorescens were not detected by the proteomic profiling analysis, the transcriptional regulation of the genes encoding for these proteins was assessed (Fig. S1b). RT-qPCR assays (Fig. S1b) showed an important transcription decrease of clsB, encoding the cardiolipin synthase B (ClsB). However, no significant expression levels difference was noted for the other studied genes upon NO 2 exposure.

Discussion
NO 2 is an important air pollutant whose effects on the human lung microbiota are poorly understood. To study the effects of atmospheric NO 2 on bacteria, the airborne P. fluorescens MFAF76a strain was exposed to either 45 ppm of gaseous NO 2 or synthetic air for 2 h, through an exposition system described previously 21 .
In this study, electronic microscopy observations revealed that exposure to 45 ppm of NO 2 affects the envelope cell morphology of MFAF76a strain (Fig. 1a). Indeed, while some control cells displayed a rough cell surface, the prevalence and intensity of these irregularities were much higher among cells exposed to NO 2 . The effect of NO 2 on the bacterial envelope is reminiscent of the study of Deupree and Schoenfisch that showed the deterioration of P. aeruginosa membranes by NO, and outlined even potential holes formation 34 , and suggests that NO 2 could alter membranes in a similar fashion. In this study, we tried to explain the emergence of the observed phenotype and to understand the effects of NO 2 on the various cell components.
Flow cytometry showed that exposure to NO 2 induces an increase in cell size (Fig. 1b,c). However, the cell size did not increase homogeneously, as the cells exposed to NO 2 displayed more high sized events (Fig. 1b,c). This increase of the number of cells of high size could be the result of an alteration of the cell division process. This alteration is however temporary as the generation time of the bacteria published in Depayras et al. were not different after exposure to NO 2 33 . The growth curve displayed a lag phase, which was at least partially explained by a significant diminution of cell viability.
To gain further insights into the observed cellular morphology modifications in MFAF76a strain, a proteome profiling approach was assessed in the absence or presence of NO 2 . Interestingly, the proteomic profiling data (Fig. 2b) revealed modifications of the amounts of several key proteins involved in the formation of the septal ring after exposure to NO 2 , supporting the results found by flow cytometry (Fig. 1b). Septal ring formation is a complex process. Briefly, FtsZ polymerizes into a z-ring at the midcell and recruits other proteins such as FtsA and ZipA that anchor the z-ring to the membrane. This complex then recruits other effector proteins that will help in the contraction of the ring and the synthesis of cell wall at the division site to prevent the rupture of the cell. Z-ring localization is regulated by the Min complex, which prevents the polymerization of FtsZ near the poles of the cell 35 . Among the proteins involved in the bacterial divisome formation, six were differentially expressed upon exposure to NO 2 . These proteins are FtsA, FtsE, FtsI, FtsZ, MraZ and MinE. The protein FtsZ can polymerize into a protofilament forming the z-ring, and allows the recruitment of other effector proteins. This polymerization is inhibited at one cell pole by the MinCD complex, which is then separated by MinE, causing MinCD to reform at the opposite pole, only to be separated by MinE again. This oscillation of MinCD and MinE restricts the polymerization of FtsZ to the center of the cell. FtsE is part of the transmembrane FtsEX complex, an ABC transporter-like protein required for cell division in low salt environment 36 . This complex is located to the septal ring through FtsE, while FtsX has multiple functions, including sequestration of FtsA to prevent its dimerization and recruitment of proteins involved in peptidoglycan hydrolysis 37 . The increase in FtsZ, MinE and FtsE quantities after NO 2 exposure suggests the activation of regulatory pathways promoting cell division. However, another essential protein, FtsA, is under-expressed by threefold under nitrosative stress. FtsA is partly responsible for the attachment of the z-ring to the cell membrane as well as the recruitment of other downstream effector proteins necessary for cell constriction. Through its ability to bind FtsZ and the membrane, FtsA plays an important role in FtsZ dynamics. FtsA acts as a membrane anchor for FtsZ, but also tends to destabilize FtsZ polymers, and an overproduction of FtsA impairs cell division. Furthermore, FtsA is also able to self-polymerize which could inhibit its functions. This process could not be totally assumed since literature reports seemingly contradictory results on this subject [37][38][39] . Therefore, several hypotheses could explain the sharp decrease of the quantity of FtsA inside the cell after exposure to NO 2 . Lowering the production of FtsA could be a way for the cell to improve the septal ring stability by promoting FtsZ polymerization. FtsA could also be degraded after failing to attach to a membrane disorganized by NO 2 . Finally, FtsI abundance is also decreased by threefold. This protein is responsible for the synthesis of peptidoglycan at the septal ring location. This biosynthesis allows the contraction of the membrane without destroying the cell envelope. Normal localization of FtsI at the septum requires FtsA among other proteins. Since most proteins involved in the septal ring division are bound to the inner membrane, either directly or through their interaction with other proteins, perturbations of the membrane integrity could modulate their efficacy.
Furthermore, several proteins involved in the synthesis of peptidoglycan are differentially produced after exposure to NO 2 (Fig. 2b). Thus, GlmS is overproduced while GlmM and GlmU are both underproduced in response to NO 2 exposure. GlmS and GlmM reaction products are not specific to the peptidoglycan pathway, so GlmU decrease is a better indicator of the effect of NO 2 regarding peptidoglycan synthesis. Peptidoglycan synthesis requires the presence of two precursors, uridine diphosphate N-acetylglucosamine (UDP-GlcNAc) and uridine diphosphate N-acetylglucosamine (UDP-Mur-NAc) 40 . Furthermore, MurC, involved in the addition of the pentapeptide moiety to UDP-Mur-NAc is also underproduced. Finally, FtsI, responsible for peptidoglycan cross-linking at the site of the division septum, and DacA, which mediates peptidoglycan cross-linking are less abundant 41,42 . These results show that peptidoglycan synthesis is underproduced upon exposure to NO 2 through the limitation of the synthesis of its precursor. Furthermore, a decrease in DacA abundance can alter bacterial morphology, by producing more irregular cells 43 . The decreased DacA quantity after exposure to NO 2 could therefore contribute to the cell surface irregularities observed by electron microscopy. The morphological alterations could also be linked to modifications of the membranes. To ascertain if the cell surface and elongated cells phenotypes are associated to membrane integrity alterations, membrane integrity measurements were evaluated in MFAF76a strain in the absence or presence of NO 2 (Fig. 3a, b). The integrity data showed that membranes of 42.5% of the cells exposed to NO 2 were completely permeant to PI, while being only weakly labelled by SYTO9™. Meanwhile, 18.2% of the cells exposed to NO 2 showed strong labelling with both PI and SYTO9™ (Fig. 3b). This co-labelling indicates that the cells where less permeant to PI, suggesting their membranes had a higher level of integrity than the cells entirely permeant to PI. These values, compared with the respective 2.5% and 7.1% of the control cells, suggest that the cell morphology alterations observed by electron microscopy could be correlated to a decrease of the membrane integrity. These results prompted us to also explore the membrane fluidity which was quantified by fluorescence anisotropy using the DPH probe (Fig. 3c). Exposure to NO 2 decreased significantly DPH fluorescence anisotropy in MFAF76a, indicating an increase of membrane disorganization. This increase could be the result of alterations to the membrane phospholipids, as well as damages to membrane proteins caused by NO 2 and its potential derivatives. Together, these results confirm that the membranes of the MFAF76a strain were strongly altered upon NO 2 exposure.
The proteomic profiling analysis of proteins involved in the synthesis of membrane components revealed that these pathways are also affected by NO 2 (Fig. 4b). UDP-GlcNAc, involved in an early step of PG synthesis, is also the precursor of lipid A, a part of lipopolysaccharide (LPS) when associated with a polysaccharide. The first stage of the LPS biosynthesis pathway is the conversion of UDP-GlcNAc into Kdo 2 -lipid A by nine enzymes 44 . The protein LpxC, which is strongly, overproduced after exposure to NO 2 , catalyses the first committed step of this conversion. Synthesis of Kdo 2 -lipid A also requires the addition of two CMP-3-deoxy-D-manno-octulosonic acid (CMP-Kdo) residues 44 . KdsC is the second of the three enzymes of the subpathway synthetizing CMP-Kdo from D-ribulose-5-phosphate 45 . This enzyme is also overproduced in cells exposed to NO 2 . Taken together, these results suggest that NO 2 had an effect on synthesis of the peptidoglycan and lipopolysaccharides precursor UDP-GlcNAc.
Furthermore, the quantity of several proteins involved in fatty acid biosynthesis (FAB) was modulated in response to NO 2 exposure (Fig. 4b). ACC is a four-subunit complex catalysing the first step of FAB in Pseudomonas 46 . Three of these subunits, AccA, AccC and AccD, were underproduced in the cells exposed to NO 2 , suggesting decreased FAB. Two proteins involved in the elongation of fatty acid, FabA and FabG, are overproduced after NO 2 exposure (Fig. 4a). Taken together, these results indicate that NO 2 exposure altered FAB pathway in MFAF76a strain. The overproduction of FabA could cause the increase of cis-unsaturated fatty acids, leading to an increase of membrane fluidity. However, the increase of membrane fluidity is unlikely to be caused by increase of cis-UFA production, since in a previous studies we showed that NO 2 exposure did not induce any significant alteration of the degree of saturation nor in cis-trans isomerization ratio of fatty acids 33 . Moreover, FadE, which catalyzes the first step of fatty acid beta-oxidation, was underproduced after exposure to NO 2 47 , indicating that the beta-oxidation of fatty acids is also impacted by NO 2 and may be downregulated. Interestingly, despite the modification observed in membrane components biosynthesis, the bacterial colonies displayed no observable differences between the control and the cells exposed to NO 2 .
The proteomic profiling study showed that gaseous NO 2 have an important effect on proteins synthesizing most of the envelope components. However, abundance variations on proteins involved in phospholipid biosynthesis was not detected with this approach. To assess whether bacteria modulate their phospholipid biosynthesis in response to NO 2 , the transcription of genes encoding for enzymes involved in this process were analyzed by RT-qPCR. No significant differences in expression levels of genes involved in phospholipid biosynthesis between control cells and NO 2 exposed cells were found except for clsB, the expression of which was reduced by eightfold. clsB encodes for the cardiolipin synthase B protein, involved in the last step of the cardiolipin biosynthesis pathway, and is only active in stationary phase of growth 48 . Despite the presence of three cardiolipin synthase genes in the MFAF76a strain, no cardiolipin was characterized in the lipidomic profile of the MFAF76a strain 33 . However, it is possible that the absence of cardiolipin could be due to the methods used to extract and analyze the cardiolipin profile. Cardiolipin is known to have an important role in membrane structure, by increasing membrane fluidity or by producing a favorable environment for the correct functioning of membrane proteins 49 . Cardiolipin is also known to be localized preferentially at the poles of the cell. Therefore it plays a role in determining the localization of proteins involved in the septal ring formation such as the Min complex by favoring its biding to the membrane at the cell poles 50,51 . These data show that exposure to NO 2 affected the cardiolipin biosynthesis pathway, which might lead to a reconfiguration of the bacterial envelope, in accordance with UGP loss upon NO 2 exposure.
Overall, these results show that NO 2 exposure of MFAF76a leads to an augmentation of the cell size, which may be linked to divisome formation impairment. Furthermore, NO 2 induced envelope alterations and decrease membrane integrity. These alterations can be partially explained by the modification of numerous components of the cell envelope, including peptidoglycan, lipid A and phospholipids. Further experiments on other bacterial strains would permit the distinction between the general effect of NO 2 on bacteria and the effects that may be specific to this strain.

Data availability
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request. Supplementary Information.